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Abstract 


The first deep field images from the James Webb Space Telescope (JWST) of the galaxy cluster SMACS J0723.3- 
7327 reveal a wealth of new lensed images at uncharted infrared wavelengths, with unprecedented depth and 
resolution. Here we securely identify 14 new sets of multiply imaged galaxies totaling 42 images, adding to the five 
sets of bright and multiply imaged galaxies already known from Hubble Space Telescope data. We find examples 
of arcs crossing critical curves, allowing detailed community follow-up, such as JWST spectroscopy for precise 
redshift determinations, and measurements of the chemical abundances and of the detailed internal gas dynamics of 
very distant, young galaxies, One such arc contains a pair of compact knots that are magnified by a factor of 
hundreds, and features a microlensed transient. We also detect an Einstein cross candidate only visible thanks to 
IWST's superb resolution. Our parametric lens model is available through the following link (https:/ /ww. 
dropbox.com/sh/gwup2Ivks0jsqe5 / AAC2RRSKceÜaX-IIFCc9vhBXa?dl—0) and will be regularly updated using 
additional spectroscopic redshifts. The model is constrained by 16 of these sets of multiply imaged galaxies, three 
of which have spectroscopic redshifts, and reproduces the multiple images to better than an rms of 0!'5, allowing 
for accurate magnification estimates of high-redshift galaxies. The intracluster light extends beyond the cluster 
exhibiting large-scale features that suggest a significant past dynamical disturbance. This work 
s a first taste of the enhanced power JWST will have for lensing-related science. 


Unified Astronomy Thesaurus concepts: Galaxy clusters (584); Strong gravitational lensing (1643); Dark matter 


atrick Kamieneski!! ®, Nicholas Fo 


(353); High-redshift galaxies (734); Galaxies (573) 


1. Introduction 


At long last, the first deep images from the James Webb Space 
Telescope (JWST) were delivered on 2022 July 11. The target is 
the central region of SMACS J0723.3-732 (SMACS 0723 
hereafter; za = 0.39), which is part of the southern extension of 
the MACS sample (Ebeling et al. 2010; Repp & Ebeling 2018). 
SMACS 0723 is a massive strong-lensing (SL) galaxy cluster 
which is rich in galaxy images distorted by the gravitational 
lensing effect, as was seen in recent Hubble Space Telescope 
(HST) imaging taken as part of the Reionization Lensing Cluster 
Survey (RELICS; Coe et al. 2019). The Planck-derived mass is 


S Original content from this work may be used under the terms 
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high 8.39 x 10" Mo, and there were nine galaxies at z > 5.5 and 
two candidates at z=7 based on their photometric redshift 
estimation (Salmon et al. 2020; Strait et al. 2021). Three far- 
infrared sources have also been detected in the field using the 
Herschel Space Observatory (Sun et al. 2022). 

A lens model for this cluster based on these previous HST 
images was published recently (Golubchik et al. 2022). This 
model uses the Light-Traces-Mass (LTM) method, which 
assumes that the galaxy light traces the underlying stellar and 
dark matter (DM) but does not assume a parametrized mass 
function for the DM (Broadhurst et al, 2005; Zitrin et al. 
2009, 2015). In Golubchik et al. (2022), the authors identify 
five strongly lensed galaxies and derive the spectroscopic 
redshifts for three of those systems using publicly available 
ESO Very Large Telescope (VLT) Multi Unit Spectroscopic 
Explorer (MUSE; Bacon et al. 2010) data. Two additional 
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Figure 1. Image stamps depicting the 19 image multiplici 


in the SMACS 0723 field, as labeled, The publicly available color imag 


‘of SMACS 0723 is depicted, 


which is valuable for the identification of image family members by their similar colors and morphological components. We also require that the model be able to 


reproduce the positions of the members of each image family. A 5” bar is shown for refe 


models are also available, one based on Lenstool (Jullo et al. 
2007; Sharon et al. 2022) and the other on Glafic (Oguri 
et al. 2013). These models, however, do not use the spectro- 
scopic information from MUSE. 

This first set of images of SMACS 0723 from JWST 
represents a milestone for not only astronomy, but science in 
general. As with its predecessor, HST, this first image from 
IWST reveals the distant universe with incredible detail. Some 
of the most prominent arcs detected with the HST RELICS 
program now show multiple unresolved substructures in JWST 
previously unseen (see for instance system 4 in Figure 1). 
These substructures facilitate the identification of systems of 
strongly lensed images, where an image system consists of all 
the gravitationally lensed counter images from the same 
background galaxy. The greater sensitivity of JWST to small 
flux fluctuations and improved spatial resolution enable 
searches for the fainter caustic or microcaustic crossing events 
of stars at z < 2 such as Icarus (Kelly et al. 2018), or Warhol 
(Chen et al. 2019; Kaurov et al. 2019), stars between 2 < z < 5 
such as Godzilla (Diego et al. 2022), stars at 26 such as 
Earendel (Welch et al. 2022) or even further up to the first stars 
as proposed in Windhorst et al. (2018). 

JWST is already providing breakthrough results in the study 
of the most distant galaxies. Pointing JWST toward gravita- 
tional lenses results in a more powerful combined telescope 
with an effective diameter a factor V[ji| times larger, where |j1| 


m 


is the absolute magnification of the lensed background object. 
For galaxies, typical magnifications can reach factors of a few 
tens, making the combination JWST+SL cluster similar to a 
space telescope analog to JWST but with a 20-30 m diameter. 
In the context of high-z object detection, the lensed source 
counts are initially lower than a blank field due to the reduction 
in search volume from magnification. At the same time, the 
blank field luminosity function flattens with higher redshift, 
allowing the number of lensed galaxies to overtake the number 
of blank field galaxies (e.g., Mahler et al. 2019; Salmon et al. 
2020; Pascale et al. 2022). 

For much smaller background objects like magnified stars 
near caustics or in general for subparsec structures near 
caustics, the underlying magnification factor can be of order 
~1000, This translates into effective apertures for JWST 
of © 200 meters! Even without the ultrahigh magnification 
boosts near the critical lines, SL has allowed us to probe several 
magnitudes deeper than blank fields, i. e., down to rest-frame 
UV luminosities Muy S —13 magnitudes (Bouwens et al 
2017; Livermore et al. 2017; Atek et al. 2018; Ishigaki et al. 
2018; Vanzella et al. 2021; Bouwens et al. 2022a, 2022b) 
and stellar masses M, Z 10°M,, (Bhatawdekar et al. 2019; 
Kikuchihara et al. 2020; Furtak et al. 2021; Strait et al. 2021). 
With the redder wavelength range, greater sensitivity, and 
resolving power of the JWST, we can expect many more 
galaxies at even higher redshifts to be detected in the coming 


THE AsTkOPHYSICAL JOURNAL LETTERS, 938:L6 (Opp), 2022 October 10 


months. In order to fully characterize them and study their 
physics, it is beneficial to have different models which are 
constructed in independent works. After this paper was 
submitted to the archive, we became aware of two new lens 
models based on Lenstool which include boh MUSE 
spectroscopic redshift constraints and JWST imaging (Caminha 
et al. 2022; Mahler et al. 2022) 

In this paper we identify 14 new image families using the 
TWST data and present a new parametric SL model. The paper 
is organized as follows: in Section 2 we detail the NIRCam 
imaging and photometry needed to provide constraints for the 
lens model, which in turn is discussed in Section 3. In Section 4 
we present and discuss the results. The conclusion appears in 
Section 5. Throughout this work we use a ACDM cosmology 
with Q4,— 0.3, 04 — 0.7, and Hy — 70km s^! Mpe~'. Unless 
otherwise stated, we use AB magnitudes (Oke & Gunn 1983), 
and errors correspond to 10. 


2. Observations and Data 


Observations with the Near Infrared Camera (NIRCam; e.g., 
Rieke et al. 2005) aboard the JWST were executed as part of 
the Directors Discretionary Time on 2022 June 6 (PI. 
Pontoppidan; Program ID 2736). Exposures were taken in 
FO90W, FI50W, and F200W in the short wavelength (SW), 
and F356W, F277W, and F444W in the long-wavelength (LW) 
channels, totaling 12.5h of integration time. Near InfraRed 
Spectrograph (NIRSpec; Jakobsen et al. 2022) observations 
were made in the same field as a part of the same program. Two 
exposures were taken in each of the FI70LP and F290LP 
filters, totaling 11787 sec of integration time in each filter. 
SMACS 0723 was also observed using the Mid InfraRed 
Instrument (MIRI; Bouchet et al. 2015; Rieke et al. 2015) and 
the Near InfraRed Imager and Slitless Spectrograph (NIRISS; 
Doyon et al. 2012). The JWST data analyzed in this work can 
be found on MAST." 

The raw data from all four instruments were processed using 
the IWST science calibration pipeline which performed the 
background subtraction, flat-fielding, correction for cosmic ray 
hits, correction for image distortions, repixelization by the 
drizzle approach onto a common astrometric reference frame, 
coaddition to make the mosaic, and combination into a color 
image. For the NIRCam data, we reduce the raw data products 
independently from UNCAL files to I2D mosaics using a custom 
pipeline based on JWST pipeline version 1.6.2 alongside the 
Calibration Reference Data System (CRDS) version 0942. The 
pipeline includes updates based on on-flight calibrations that 
account for zero-point offsets between the filters and the 
NIRCam modules (Rigby et al. 2022), critical for accurate 
photometry and derived properties. Furthermore, this pipeline 
includes additional steps over the default JWST pipeline to 
improve on its astrometry and background subtraction. For a 
detailed description of the additional steps, see Adams et al, 
(2022). 

Using the NIRCam images based on our independent 
reduction with improved calibration, we extract photometry 
for our targets of interest using forced photometry of small, 
0/32 diameter circular apertures centered on the peak flux of 
the target in each band. The photometry is aperture corrected 
for a point-spread function using the simulated point-spread 
functions (PSFs) from WebbPSF (Perrin et al, 2014). This 
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Table 1 
List of Are Systems 
1D RA Decl pom 
120000 120000 
m 0723218012 732703467 145 
12 3430 -73:27:17:17077 — 145 
13 PLATA -73:27:31.329 — qa 
21 213288 —73:27:03369 138" 
22 21838 —73:27:18.337 138" 
23 208 —73:27:31.378 — L3 
sr 07:23:19.3406 — —73:26:54.554 181 
11.65-1.91] 
ar 07:23:19.6709 — 132718562 ‘ 
ix 0723080861 — —732734507 — 66500 , 
3a 0723:7.6080 — —732717.228 , 
m 072313007 — —132716313 — 664 0 202 
[1.94-2.48] 
42 Bo —73:27:30.116 — 24653 ^ 
43 152214 —732655308 — 22015 D 
54 17715 -73:27:06.454 — 14$ 
52 174012 — -732109402 — 4$ 
53 170777 —73:27:36422 — 0 
5A 175037 — —T32708213 — qm 
5B 176067 -73:27:08.073 — 14$ 
61 206382 — —73:27:06.170 — 26610 m 
1127-202] 
62 208612 -73:27:10.805 — 13T'12 ^ 
63 202857 — —73:27:39.168 , 
7i 109765 — —73:26:54.832 627 
1295-838] 
72 07:23:10.7520 — 13:26:56,676 , 
73 07:23:11.9016  —73:26:49389 3.66208 D 
tm 07:23:12.6279 — —73:27:36.552 593 
(24.67) 
82 0723:1.9309 — ~73:27:19.101 ^ 
83 072345.38 — —73:2648.223 — 6260 s 
on 0733027404 73:27:17.641 — 20994 290 
12.57-3.14] 
92 132719 —73:27:32.036 , 
93 15.0263 — ~73:26:53.104 D 
10.1" 17.6727 ~73:27:06.073 
10.2" 172165 — —73:27:09.979 
103" 169616 — -732736214 — 1190 
m 144676 — —73:27:21.001 1.66 
11.22-198] 
ua 0723046783  —73:27:25.882 " 
121 07:23:17.3278 -73:26:56.724 1.56 
11.52-1.71] 
07:23:15.5468  —73:27:15.540 — 675005 i 
072362084 — -132733244 — 67908) , 
134 07:23:19.1576 — —732656007 — 173133 234 
1229-277] 
173344 — —73:27:14.825 ' 
176167 -73:27:41.842 — 2.69508 D 
19.8353 — —73:27:15.697 , 
141! — 072312422 — —T32729351 183 
[1.68-2.03] 
142" 0723:24965 — -73:27:31.974 ^ 
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Table 1 
(Continued) 
D RA. Dec p m 
320000 3200.0 

3 0723125821 — 7E 7 

15419 0723167116 — —73:26:55.135 

152 072344745 — —73:27:16.356 6 597288 

153 0723154182 — ~73:27:32.208 2.300008 , 

161 072369067 -73:27:09.328 — 02909 LoT 
1082-104] 

162 0723169206 — ~73:27:10.256 ^ 

74 072347.975 —73:27:20.571 — 2450 172 
[1.68-1.92] 

172 07:23:17.7931 — 132726849 a 

173 072319770 — 132650524 , 

18.1 0723472313 — —73:27-02.060 

182 07:23:16.040 — 732713482 

183 0723163581 — ~73:27:32.183 ^ 

GLI' — 072321876 -73:27:41.814 Lal 
[1.15-1.58] 

GL2! — 0723224453 -73:27:40.964 1.498888 ^ 

GLI” — 072322478 — —732741727 — 1490 £ 


Einstein — 0723027975 — -732708814 


Notes. Column 1: ID; Columns 2 & 3: R.A. and decl.; Column 4: Redshift, For 
systems 1, 2, and 5 we quote the spectroscopic redshift from MUSE 
(Golubchik et al. 2022), and note that the NIRISS spectrum for system 1 
confirms the stated redshift. Column 5: The redshift of the system as predicted 


from the SL model. Candidate images whose identification is not secure are 
marked with “2.” Note that these were therefore not used in the minimization 
? The fist five image systems are drawn from Golubchik et al. (2022). 

P These are the photometric redshift estimates using the six bands of NIRCam 
imaging, 

© These redshifts are spectroscopic values from Golubchik et al, (2 
“ These systems were not used to constrain the parametric SL model, but were 
confirmed in the WSLAP+ model. 

© Mahler et al. (2022) find a spectroscopic redshift of 1.99 for this system using 
MUSE, but this is not included in the model. 


2. 


information is then fed into the photometric redshift code 
LePhare (Arnouts et al. 1999; Ilbert et al. 2006). LePhare 
operates by fitting a large grid of galaxy spectral energy 
distributions (SEDs) to the photometry provided. We run 
LePhare with a suite of galaxy templates from (Bruzual & 
Charlot 2003), allowing for E(B-V) values between 0 and 1.5 
(Calzetti et al. 2000), redshifts between 0 and 12, and applying 
the (Madau 1995) treatment for absorption from the IGM. The 
redshifts provided in Table 1 are the redshift of peak 
probability and errors generated from a x^ grid produced by 
LePhare based on the templates, redshifts, and dust 
attenuation combinations available. Several entries are missing 
values as a result of various sources of contamination or by a 
lack of detection. We refer to the Appendix for more details. 
We note that these NIRCam data can be compared with bluer 
data using HST taken as a part of the RELICS program (PI: D. 
Coe), and the Spitzer Space Telescope (PI: M. Bradac). 
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Figure 2. Close-up view of system 5 with the z = 1.43 critical curve marked in 
the dashed line on this FO90W--FIS0W--F200W composite color image. Two 
compact sources in the inset, knots A and B, have a separation of ~0/"16, and 
have sizes that are consistent with the width of the PSF. The knots bracket the 
critical curve that must pass in between them at an angular separation 
of ~ 008, We estimate high magnification factors as a result of this close 
proximity to the critical curve of |j1| ~ 750 for each source. A red knot is also 
detected at a position slightly offset from knot A, Tr, which may be a 
microlensed transient. We refer to Section 4 for more information regarding 
this potential microlensed transient. 


We also made use of the available NIRSpec data set. In 
particular, we analyzed those spectra that were acquired at the 
positions of our image systems to confirm our image system 
designations. We refer to Section 3.1 for more details. We also 
consulted the MIRI images to test for achromaticity within each 
of our image systems, for those systems that were detected. 

Overall, the arcs in this field make up a rich tapestry of 
distorted images typical of lensed sources and at a level of 
detail that has never been seen before. The increased spatial 
resolution reveals lensing constraints even down to the small 
substructures within the arcs, since many of these small-scale 
sources are equally multiply imaged, such as the double are in 
system 5 (Figures | and 2). Meanwhile on the scale of the 
cluster, thanks to the remarkably dark infrared sky background 
of JWST images, we can get a high signal-to-noise view of the 
intracluster light (ICL). This cluster has an ICL that is 
elongated and nonuniform, We refer to Section 4 for a 
discussion of how the ICL's relatively unique features suggest 
that this cluster is not dynamically relaxed. 


3. SL Modeling of SMACS 0723 


We generate a new SL model exploiting the JWST data in 
which we identify the new sets of multiple images, as described 
in Section 3.1, Cluster member galaxies used in the modeling 
are described briefly in Section 3.2. The modeling method is 
described in Section 


n 
3. 


31. 


Yew Multiple Images in SMACS 0723 


We identify new lensing constraints by taking advantage of 
the superior quality of the JWST images. We start by adopting 
the first five multiple image systems published in Golubchik 
et al. (2022), following their numbering scheme. These systems 
serve also as a model-independent guide to search for new 
systems, For example, galaxies at close projected separations 
and at similar redshifts typically maintain similar s 
distributions when distorted by a gravitational lens. Moreover, 
among the galaxy images, counter images with negative and 
positive parities are in general symmetric to one another with 
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respect to the critical curve. After visual inspection of the new 
IWST/NIRCam image, and corroboration by available photo- 
metric and/or spectroscopy redshift constraints and model- 
predicted locations, we identify 14 new sets of multiply lensed 
galaxies which are reported in Table 1. 

Fortunately, there are spectroscopic redshifts available for 
some of the image systems. The redshifts of image systems 1, 
2, and 5 are z= 1450 -- 0.001, z= 1.378 40.001, and z= 
1425 + 0.001, respectively and were measured with MUSE 
(Golubchik et al. 2022). Note that while JWST/NIRSpec 
spectroscopic redshifts of 35 galaxies have been reported in this 
field with redshifts as high as 8.3 as reported by the JWST 
press release, an investigation of the NIRSpec spectroscopy 
uncovered spectra for only one of our 42 images from our 19 
image families: image 4.1. However, an inspection of the 
spectrum in each of the spectral bandpasses did not produce 
any significant features that were present in multiple visits and 
therefore did not yield a redshift. While valuable, further 
spectroscopic redshifts are however not crucially necessary for 
securing our image system identifications which rely on a large 
set of image multiplicities vetted in several ways. We note that 
NIRISS spectroscopy of SMACS 0723 is also available but not 
used in this study. Its redshift constraints will inform future lens 
models. 


3.2. Cluster Member Galaxies 


For our model we use the same cluster member galaxy 
selection as in Golubchik et al. (2022), which was based on the 
red sequence shown therein (their Figure 4). In addition, we 
add here three to four galaxies not used in the initial model 
by Golubchik et al. (2022) which lie near the bright star 20" 
north of the BCG (which hindered somewhat their previous 
detection) and are included as they may affect the reproduction 
of nearby counter images. A revision of this selection using 
spectroscopic data (such as e.g., the NIRISS data), is deferred 
to future work. 


3. SL Modeling Method 


The LTM model presented in Golubchik et al. (2022) 
required a very strong external shear, indicating that the 
cluster's central mass distribution is highly elongated. Since the 
LTM approach is limited in the intrinsic matter ellipticity it ean 
immodate, and for comparison, we make use here of a new, 
fast parametric method that we recently constructed (A. Zitrin; 
2022, in preparation; the method can also be referred to as 
analytic, i. e., it is not limited to a grid’s resolution). The 
method is similar in nature to other parametric lens modeling 
techniques such as Lenstool (Jullo et al. 2007), Glafic 
(Oguri et al. 2010) or GLEE (Halkola et al. 2006; Grillo et al. 
2015). In such parametric methods, dark matter halos can be 
elliptical which more easily accounts for the required elongation 
(which is also evident from the distribution of arcs). 

The models in our new parametric method have two main 
components, First, the cluster members galaxies are modeled with 
double pseudo elliptical mass distributions (PIE; Elfasdóttir et al. 
2007) which are assumed spherically symmetric with the exception 
of the BCG, and are defined following the prescriptions of Jullo 
et al. (2007) and Zitrin et al. (2013): 


Pascale et al. 


humum) o. oy 


Where the ø is the velocity dispersion, re is the cutoff radius, 
Teore ds the core radius, and L' is the typical luminosity of a 
galaxy at the cluster redshift. We fix 7, = 0.2 kpe for all 
dPIEs, use an L” equivalent to a galaxy of misi4y = 22.13, and 
assume a constant mass-to-light ratio (a =0.5), while o” and 
rh are left free to be optimized by the model. 

The second modeling component is the set of cluster DM 
halos. These can in principle be modeled either as elliptical 
Navarro-Frenk-White (NFW) profiles (Navarro et al. 1996), or 
as pseudo elliptical mass distributions (PIEMD; e.g., Monna 
et al. 2015); or dPIEs. The third and last component which can 
be added if necessary is a two-component external shear. 

The method is similar to our previous parametric imple- 
mentation outlined in Zitrin et al. (2015) which has been well 
vetted, but is not limited to the assigned grid resolution. The 
new improved version used here has been already applied to 
various clusters and has also been tested on simulated clusters, 
accommodating both image- and source-plane minimization. 
The minimization of the model is done via a Markov Chain 
Monte Carlo (MCMC) with a Metropolis-Hastings algorithm 
(e.g, Hastings 1970). We include annealing in this procedure, 
and the chain typically runs for several dozen thousand steps 
after the burn-in stage. Errors are calculated from the same 
MCMC chain. 


34. SMACS 0723 SL Model 


For modeling SMACS 0723 by this method, we constrain 
the model using 16 of the 19 identified image systems and a 
total of 48 images (two sets of multiple images are used from. 
system 5; see Table 1), including spectroscopic redshifts for 
three systems and leaving all other system redshifts to be fit by 
the model. We represent cluster members as dPIE profiles, as 
explained above, and the cluster DM halos as PIEMDs (i.e., the 
main free parameters for each halo aside from the ellipticity and 
its position, are a core radius and a velocity dispersion). We use 
two DM halos, centered on the two central galaxies 
respectively (marked by “o” symbols in Figure 3). We do not 
incorporate an external shear. Since galaxies can deviate from 
the assumed scaling relations, we leave the four brightest 
galaxies to be freely weighted (marked by rings in Figure 3), 
which means each of their velocity dispersions c is scaled by a 
weighting parameter optimized by the model (see Table 2; G1- 
G4). The BCG's ellipticity and position angle are also left free 
to be optimized by the model. In total, there are 29 free 
parameters in the model: the normalization for the galaxy 
velocity dispersion o", the normalization for the galaxy cutoff 
radius rå, four free parameters for each of the two dark matter 
halo PIEMDs, four galaxy weights, the BCG ellipticity and 
position angle, and the 13 redshifts for systems without 
spectroscopic redshifts. For the lensed sources, each multiply 
imaged system should converge to a single position in the 
source plane. To enforce this constraint, we minimize the 
differences between the image positions in the source plane 
following the prescription of Keeton (2010), which usually 
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Figure 3. The critical curve at z = 1.45 from the SMACS 0723 lens model is overlaid as a color composite image generated from the six existing JWST NIRCam 
filters depicting the central 1.5" by 1/ of the cluster. The model incorporates strong-lensing constraints from the 19 image families identified in total, 14 of which are 
discovered by this work using the JWST data). The two “+” symbols mark the positions of the two dark matter halos used in this model, and the four magenta rings 
mark the galaxies whose weight is lef free in the model. Image 2.2 is rare for containing > 10 bright knots located at large projected galactocentric radii (inset). 


Table 2 
SL Model Optimized Parameters 

E 35 Tus Tok) 
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Notes. Column 1: Object; Column 2: Ellipicity: Column 3: Position Angle: Column 4: Velocity Dispersion. Column 4: Cutoff Radius. All errors are given as the 
68.3% Confidence Interval 
* The scaling relations for all galaxies 


P Weights for the four galaxies allowed to deviate from the scali 


converges to a solution relatively quickly and without loss of 
quality, as shown by Keeton (2010) in their successful Hubble 
Frontier Fields (HFF; Lotz et al. 2017) lens models, Given the 
48 lensed images, this constraint amounts to 62 SL constraints 
on the model (two for each lensed image with the exception of 
one image per multiplicity), which gives a total of 33 degrees 
of freedom for the model 


4. Results and 


ussion 


The best-fit SL model, which is the one for which y? 
is minimized, is presented in Figure 3. Critical lines are 
overlaid for z, = 1.45 (i. e., the redshift of image system 1). The 
figure illustrates the elongated configuration of SMACS 0723. 


relations, as labeled in 


igure 3. 


The model was minimized using a positional uncertainty of 
Apos = 075, resulting in a x? of 159.4 (x? = 4.8) in the image 
plane, and an rms of 0/48 in reproducing the positions of the 
multiple images; the optimized model parameters can be found 
in Table 2. 

The effective Einstein radius at a given redshift is computed 
as @,=JA/z, with A being the area enclosed within the 
critical curves. At z=2, Øp=18”4+1”8, and the mass 
contained inside the critical curve is (5.91 +0.83)x 10^ Ma, 
which is overall consistent with the values of 0; — 16"9 + 2”, 
nd M = (4.15 + 0.58) x 10" Mo obtained in the HST LTM 
model (Golubchik et al. 2022). The uncertainties on the 
Einstein radius and mass are similar for each of the models and 
are limited by the systematics. 
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The JWST/NIRCam images uncover fine details of the faint 
substructures within the ares. In particular, two compact 
sources in system 5, knots A and B (5A and 5.B in 
Table 1), are spatially resolved on opposite sides of the critical 
curve, constraining the position of the critical curve in our lens 
model (Figure 2). The image pair has a close angular separation to 
the critical curve of 0"08, implying extremely high magnifications 
for each image and making it an ideal system in the search for 
caustic crossing events, Near to the critical curve, magnification 
follows an inverse relationship with distance ji = j1o/Deriy Where 
Ho is dependent on the slope of the lensing potential and can be 
measured from the lens model (Diego 2019; Welch et al. 2022). 
Uncertainty in jig is accounted for by sampling the model MCMC 
chain, while the uncertainty in Deg is assumed to be the 
uncertainty in the centroid position of the images. Our parametric 
model implies a magnification of jj = 737* 137^ for each image, 
whose positions are reproduced to within 0702 by the model. 

Magnification near the critical curve is highly sensitive to 
systematics, and may vary significantly between modeling 
approaches (Meneghetti et al. 2017). Meneghetti et al. (2017) 
demonstrated that the uncertainty in magnification could reach 
greater than 30% at ji» 10 across all modeling approaches. 
Hence statistical errors in magnification from the model are most. 
likely to underestimate the true error which is dominated by 
systematics. Assuming these images are unresolved, this magni- 
fication could imply parsec or even subparsec sizes for this object. 
Given the are’s close proximity to the cluster center, however, 
microlensing from intracluster stars could work to smooth out the 
critical curve, placing an upper limit on the persistent magnifica- 
tion possible (Venumadhav et al. 2017; Dai 2021). 

Consider again the image pair consisting of knots A and 
B in Figure 2. There is a red knot offset from knot A that is 
detected only in some of the NIRCam bands, called Tr (see 
Appendix B). We present the argument that Tr is neither knot 
A nor knot B. Assuming knots A and B are counter images, the 
pair will always be detected together in any given band. Since 
the separation of these mirrored knots is about 0.16", which is 
larger than the FWHM of any JWST filter in this data set 
(F444W, FWHM = 0.145", is the largest), JWST should be 
able to spatially resolve this pair. We detect knots A and B in 
all the SW filters, but detect only a single peak in each of the 
LW filters near to the position of knot A (see Appendix B and 
Figure 6). Since the image pair is not detected together through 
the LW filters, this single peak must refer to Tr. Hence, Tr is 
distinct from knots A and B. The lack of a counterimage, and 
its placement in close proximity to the critical curve, suggest Tr 
is a microlensing event. However we cannot rule out that it is a 
red foreground galaxy. Given the very red colors and small 
angular sizes involved, confirmation of Tr as a transient would 
require JWST follow-up observations, We note that there are 
no other similar dual-color detections in any other objects in 
system 5 that would suggest an image misalignment. 

The giant arcs also present some rather unusual features, 
such as system 2 with a spectroscopic redshift of z = 1.38 first 
measured by Golubchik et al. (2022). Image 2.2 has >10 
prominent clumps all of a similar red color and extending to 
high galactocentric radii (see Figure 3, inset). While the 
significant subkiloparsec structure of star formation regions is 
not unusual for star-forming galaxies at z— 1-3 (or even 
smaller than 100 pc; e.g., Johnson et al. 2017), the very large 
separations cast some doubt that the extended knots are 
identical in nature to those in what might be the plane of a 
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galactic disk. These clumps also appear to be different (redder) 
than those typically detected in clumpy galaxies (e.g., Shibuya 
et al. 2016). The physical characteristics are more consistent 
with numerical predictions for globular clusters, For example, 
Sameie et al. (2022) found in hydrodynamic simulations that 
the longest-lived globular clusters likely formed at high 
redshifts beyond the half-light radius (up to hundreds of 
parsecs away from their host galaxy). Pozzetti et al. (2019) go 
on to predict the color evolution for high-redshift globular 
clusters, but also point out that NIRCam colors are most likely 
insensitive to the formation redshift. Vanzella et al. (2017) have 
reported detections of candidate young (<10 Myr) globular 
clusters which are blue and at high redshifts of z = 3-6. Based 
on this information, it is tempting to speculate that the clumps 
seen in system 2 may be redder and older (~Gyr) counterparts 
of those young candidate globular clusters detected at the much 
higher redshifts, Based on SED fitting from 0.4 to 4.4 jum, 
Mowla et al. (2022) found evidence that the clumps were 
consistent with evolved globular clusters with ages of 
3.9-4.1 Gyr, implying formation a mere ~0.5 Gyr after the 
Big Bang. We note that another possible explanation, is that the 
red clumps fitting the description of those seen in system 2 may 
arise as a result of ram pressure stripping of dwarf galaxies in 
the dense environments near in projection to its disk (e.g. 
Mayer et al. 2006; Boselli et al. 2022). Further follow-up, such 
as detailed source plane reconstruction and photometric SED 
modeling, is needed to uncover the nature of these sources. 

A striking feature of this new JWST/NIRCam image is the 
morphology of the baryonic component responsible for the 
intracluster light (ICL), which exhibits a loop-like feature in the 
northwest component and a large lobe in the southeast 
component, We perform our own reduction of the raw images 
to mitigate the variations in sky background across chips as 
described in Section 2, although some background noise can 
also be seen in each of the four comers of the image due to 
unoptimized sky subtraction in the JWST reduction pipeline 
(Figure 4), The ICL is due to stars stripped away from their 
galaxies but still gravitationally bound to the cluster. Earlier 
work has suggested that the ICL is a good tracer of the DM 
distribution (or vice-versa) since both stars in the ICL and DM 
are expected to behave as collisionless particles and thus 
respond only to gravity (Montes & Trujillo 2018). Our initial 
SL model exhibits a similar macrostructure to the ICL as seen 
in Figure 4 and in Figure 5. The ICL also extends beyond the 
current set of SL constraints, but in combination with weak 
lensing derived from the JWST images it should be possible to 
produce a lens model that covers the entire range of the ICL in 
a way that gives insights as to the intriguing apparent 
correlation between the ICL and the DM distribution, In 
particular, the presence of these large-scale features in the ICL 
may motivate detailed simulations to examine the merging 
scenarios that constrain the time since the major merger and 
estimate the relative velocities involved. 

We note that after this paper was submitted to the archive, 
and prior to submission to ApJ, we became aware of two other 
papers on SL modeling in this field (Caminha et al. 2022; 
Mahler et al. 2022). There are some differences in the results 
presented in these other papers, but we have checked and found 
that including their results in our lens model does not alter the 
results of this study. 


THE ASTROPHYSICAL JOURNAL LETTERS, 938:L6 (Opp), 2022 October 10 


Figure 4. The intracluster light for SMACS 0723 in a F277W-+F356W 
+F444W color composite image. The ICL structure is especially clear in 
F277W and F356W, and is much less prominent in the SW bands. The 
pronounced ICL extends along the long axis of the cluster similar to the mass 
distribution from the lens model, as shown by the overlaid contours of 
projected mass density «. Typically a smoothly varying feature, in SMACS 
0723 the ICL has significant structures, including a large “loop” feature in the 
northwest component, and a large lobe-like feature in the southeast component. 
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Figure 5. The projected mass density of SMACS 0723. We show «, the surface. 
mass density distribution in units of the critical density for strong lensing, for 
the source redshift of system 1, at z = 1.45, plotted on a grid of constant R.A. 
and decl (dashed lines) The resulting mass distribution exhibits a clear 
ellipticity along the same axis as the ICL as shown in Figure 4. 


110.8250" 


5. Summary 


We presented a new JWST parametric model for the massive 
cluster SMACS 0723, incorporating new multiple image family 
constraints identified in the JWST press release color image and in 
the NIRCam data set. Our model builds upon the existing lens 
models (Golubchik et al. 2022) by increasing the SL constraints 
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with 14 newly identified image families, which represents a factor 
of three improvements over the HST model. This parametric model 
has a large impact because it provides a starting point for modeling 
SL clusters in the era of JWST. The mass map reveals a somewhat 
complex and extended mass profile and, in the detection of the 
image pair knots A and B opposite the critical curve in system 5, 
also opens a route to the study of caustic transients in the discovery 
of an apparent microlensed transient adjacent to knot A. The ICL 
has a shape that is broadly similar to that of the galaxies and the 
model, but with rare larger-scale features that appear to retain the 
memory of its dynamical history, a question which can be 
investigated further with numerical simulations. In the future, 
spectroscopic and photometric redshifts will bolster the reliability of 
our lens model, allowing also a more precise placement of the 
critical curve needed to engage in caustic transient studies. 
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Appendix A 
Are Systems in SMACS 0723 


We present the list of all multiply imaged systems below in 
Table 1. The image families are vetted in multiple ways: by 
their similar morphological components, similar colors, similar 
redshifts (photometric and/or spectroscopic), and/or consis- 
tency with the lens model, The IDs are given in the first 
column, where the “?” indicates an image member candidate, 
by which we mean that it fails one or more of the above criteria, 
or that there is more than one candidate that fits these criteria. 
The photometric redshift estimates, zyizcam are computed 
using the SED fitting code LePhare, as discussed in the main 
text, with matched-aperture photometry of the six bands of 
NIRCam imaging as described in Adams et al. (2022). An entry 
has no zwmcas Value if the lensed source was significantly 
contaminated by a stellar diffraction spike, a bright cluster 
member galaxy or other projection effects, or if it was too faint 
to be detected in two or more of the NIRCam bands. The 
model-predicted redshifts, Zmodeh and ZNiRCam are in relative 
agreement with the exception of five systems. The discrepan- 
cies in systems 3, 4, 12, and 15 can be attributed to a 
degeneracy in the SED fitting for z=6-7 and z= 1.5-2.5, 


F150W. 


F277W. 


F356W. 
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which comes as a result of the Lyman Break and Balmer Break 
both falling between F090W and FISOW at these redshifts 
respectively. The inclusion of FIISW could break this 
degeneracy in most cases, and may be necessary for any future 
surveys of the redshift distribution of background galaxies. The 
other outlier is image 16.1, for which zxicam — Zmodet > 30. 
This system is faint and elongated, and is most likely not well 
represented by the 0/32 circular aperture imposed on these 
data, in addition to being contaminated by the ICL. We also 
note that the lens models are limited by the arc systematic: 
Lens predictions may change as image system constraints 
become improve, for example by the measurement of new 
spectroscopic redshifts, The optimized paramters for the SL 
model can be found in Table 2. 


Appendix B 
Image Pair in System 5 


Knots A and B of system 5 are counter images that are 
situated opposite the critical curve and thus likely highly 
magnified. As shown in Figure 6, the image pair appears in 
each of the SW filters and is spatially resolved, with an angular 
separation of 0.16". In SW filter F200W, the microlensed 
transient, Tr, is detected as well; it is bright and offset from the 
position of knot A, resulting in the red knot seen in Figure 2. In 
the LW filters, the knots A and B drop out, leaving only the 
redder Tr visible as a single peak. 


F200W 


F444W. 


Figure 6. The image pair A and B of system 5 across all six available JWST filters. In each filter we center on the positions of A and B with rings of radius 0.16", 
which is roughly equal to their angular separation, The FO90W, FISOW, F200W filters exhibit two separate knots consistent with the marked positions of A and B. 
F200W. however, also shows the bright source Tr that is slightly offset from the position of knot A. Only a single peak is detected in each of the LW filters, which is 
also slightly offset from the position of knot A. The insets depict the F200W and F444W images with the diffuse arc background subtracted off and the position of 
knot B indicated by a “+” sign. Since knot A and knot B are counter images of one another and must appear together, the Tr object seen in F200W and the LW filters 


is instead understood to be a candidate microlensing event. 
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